Evolution of fragmented states 
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Abstract. We consider the problem of evolution of the many-body state of a weakly 
interacting system of bosons in an initially fragmented (Fock) state. We show that the state at 
5-H , any time can be expressed as a continuous superposition of an infinite number of Gross-Pitaevskii 

states. 
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Quantum mechanics incorporates in a single elegant formalism features that resemble wave- 
like behavior and particle- like behavior. The many-body field operator ip{x) for example, not 
only describes the annihilation/creation of a particle at position x, but it comes as well with a 
phase, as do classical fields describing wavelike phenomena. The phenomenology associates with 
the quantum- mechanical phase is vast: narrowing down our focus to the literature related to 
superfluidity and superconductivity, we find that phases play a crucial role in the description of 
the Josephson effect. Here we concentrate on the case of atomic Bose-Einstein condensates, in 
which Josephson effects are by now theoretically and experimentally on solid ground pQ. In this 
case, the system consists of two condensates, and the relative phase between them is assumed to 
be established previously, as in the case of metallic superconductors separated by an insulator, 
in which electrons can tunnel between the two sides before and during the cooling process. With 
atomic Bose-Einstein condensates a new type of experiment [2] is possible, if two such gases are 
connected by a weak link or put to interfere only after they have been cooled separately below 
the Bose-Einstein transition temperature. Now the initial state is a fragmented (Fock) state 
| N/2,N/2)t=o, with no well-defined relative phase between the two components. In fact, one 
of the first experiments [3] with condensed atomic gases has been along these directions, and a 
many others followed [3]. In the non-interacting case, the theory describing the measurement- 
induced phase coherence is now well-established [5]; also, information about correlations in 
low-dimensional quasicondensates can be obtained by this interference technique [6]. 

But describing interfering interacting gases which are prepared independently remains a 
challenging problem. We analyzed recently this situation in |7J. In this paper, we give more 
details about the derivation of an approximate simple expression for the many-body wavefunction 
for the evolving system. We will use the subscripts L and R to denote condensates prepared 
initially in a left- and a right- well, separated by a high enough potential barrier so that tunneling 
is not possible before an initial moment t = 0, when the condensates are allowed to recombine. 
The field operator is ip(x), x is the coordinate, and the initial state is a fragmented (Fock state), 
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where a,L,R{t = 0) = / dxQ* L n(x,t = 0)tp(x). At t = 0, the trapping and the barrier potential is 
switched off; the Hamiltonian for t > becomes (with g = 4irh 2 a/m, a > 0) 



if 



2m 



(2) 



What is the structure of the many-body state at t > 0? It is easy to convince oneself that this 
Hamiltonian does not preserve the fragmented structure Eq. (JTJ) of the initial state. Another 
approach is then needed. We expand the fragmented state Eq. ([T]) into phase states [5] 
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where c = 2 N / 2 (N/2)\/y/W. 0. 

With the notation | ^^(t = 0))n 
the phase state, we have 
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This form is not invariant either under the action of the evolution Hamiltonian. However, we 
show that within the same degree of approximation as the standard time-dependent Gross- 
Pitaevskii equation, the state at a later moment t > can be written as 
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So although a fragmented state does not have a single order parameter, as we will show below, 
it can be approximately described by an infinity of order parameters y/NQ^fa, t). This comes 
from the fact that the evolution Hamiltonian, in the limit N — > 00, has approximately zero 
matrix elements between states | ^^(t))^ with different phases tp; therefore each of these states 
is evolved independently. We start our proof by noticing that, as an immediate application of 
the Cauchy-Bunyakovski-Schwarz inequality in L 2 Banach spaces, | (<Ev(£) | 3><n(i)) | reaches a 
maximum value of 1 for tp' = tp. The quantity |jy ($„,/(i) | 3><7>(i))jv | will then be very close to 
zero for tp' ^ tp, and strongly peaked to 1 if <p' = tp, so we can write 
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where a(t) = \(d 2 / dp' 2 ) | (3y | $ v )(t) 

We now introduce [8J, at any time t and for a given phase tp, a projection operator R v (t) 
orthogonal to the subspace spanned by the ket | <£^(i)), R v (t) = 1— | ^(^(^(i) |. Then, at 



any time t, we split the field operator into a condensate part of the order of v ./Vxwavefunction 
and smaller phonon part of the order of lxwavef unction, ip(x) = ^,p{x,t) + x v (x,i), where 
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Xtp(x,t) = J dx'R, iD (x,x';t)^ ip (x'). 
It is straightforward to check that x^{x,t) satisfies the commutation relations: 



Xip{x,t),xUx' ,t) =R tp {x,x',t) 
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In the Heisenberg picture, the evolution of the phonon operator x (x,t) (we put an extra 
superscript (H) for operators written in the Heisenberg representation) is given both by the 
Heisenberg commutator and the intrinsic evolution due to the change of the modes on which 
condensation occurs: 
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We proceed now to the next step, which is to linearize the Hamiltonian ([2]) around & v $L 
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then calculate the commutator, and, keeping only the higher-order contribution (in yN), we 
g6t d 

ih—Xtp(x,t)=i% J dxR !fi (x, x'\ t)tp(x') ~ ifi J dx'R !fi (x, x' , t)$ ¥ ,(x', t). (10) 

But f Rip{x, x' , t)^ v {x\ t) = 0, therefore taking a time-derivative we obtain for the intrinsic 
time-dependence the formula 
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For the usual Heisenberg commutator of @ we use the commutation relations © and the 
expansion form of the Hamiltonian and obtain 
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As a result, we finally get for ([9]): 
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The consistency condition for the expansion tf){x) = § v (x,t) + x<p(x,t) is 
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or in other words the quantum depletion of a state (p should not acquire contributions of the 
order of the corresponding order parameter during evolution: 
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Thus, to obtain the approximate evolution for the fragmented states which satisfies Eq. (I13D we 
have to impose a time-dependent Gross-Pitaevskii evolution for each wavefunction $p(x,t), 
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ih—<S> v (x', t) = "^V 2 ^(x, t)+gN\ %{x, t) | 2 %(x, t). (14) 

It is also easy to convince oneself, that, by the same procedure, one can show that the 
approximation works for any many-body operator constructed from ip. Indeed, in the second 
quantization any number-conserving operator can be written as O = J2k ®kW \ V'L where 
Ofc [V^ , -0] is a functional in which each operator ffl and -0 enters k times. Then we use 
tp(x) = $tp(x, t) + x v (x, t), we expand in powers of x<p(x, t) and, by following a similar procedure 
as above, we calculate ih(d/dt)t(&\Okffi , V>]|\E'}t and show that the first-order contribution in 
X(p(x,t) vanishes, if Eq. (I14|) is satisfied. Therefore, in the thermodynamical limit and given 
that k *C N, we obtain 

t(*\6\*) t c^^Y,^ J dpa v (t)O k [%(t),%(t)]. (15) 
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For example, if we consider the kinetic energy operator — (h 2 /2m) J dxip^ (x)V 2 ^(x), we find 
in the same way, by expanding the field operator, that 

t(*| dx^(x)V 2 i>(xm t ~ - C ° 8 ^ 2 I J dxd^itW^x^V^^x.t). (16) 

In conclusion, we show that the evolution of a fragmented (Fock) state of two weakly 
interacting Bose-Einstein condensed atomic gases can be approximated by a many-body state 
consisting of an infinite superposition of Gross-Pitaevskii states with diferent initial conditions, 
corresponding to different initial phase differences between the condensates. 
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